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ABSTRACT 

The Sloan Digital Sky Survey (SDSS) source J102347.6-I-003841 was recently revealed to be a 
binary 1.69 millisecond radio pulsar with a 4.75 hr orbital period and a ~0.2 Mq companion. Here 
we analyze the SDSS spectrum of the source in detail. The spectrum was taken on 2001 February 1, 
when the source was in a bright state and showed broad, double-peaked hydrogen and helium lines — 
dramatically different from the G-type absorption spectrum seen fr'om 2002 May onward. The lines 
are consistent with emission from a disk around the compact primary. We derive properties of the 
disk by fitting the SDSS continuum with a simple disk model, and find a temperature range of 2000- 
34000 K from the outer to inner edge of the disk. The disk inner and outer radii were approximately 
10^ and 5.7x10"'^° cm, respectively. These results further emphasize the unique feature of the source: 
it is a system likely at the end of its transition from an X-ray binary to a recycled radio pulsar. The 
disk mass is estimated to have been ^ 10^^ g, most of which would have been lost due to pulsar wind 
ablation (or due to the propeller effect if the disk had extended inside the light cylinder of the pulsar) 
before the final disk disruption event. The system could undergo repeated episodes of disk formation. 
Close monitoring of the source is needed to catch the system in its bright state again, so that this 
unusual example of a pulsar-disk interaction can be studied in much finer detail. 
Subject headings: binaries: close — stars: individual (J102347. 6-1-003841) — stars: low-mass — stars: 
neutron 



1. INTRODUCTION 

The Sloan Digital Sky Survey (SDSS) source 
J102347.6-h003841 (hereafter J1023) is located well off 
the Galactic plane at 6 = 45.8°. kt V = 17.5 and Kg = 
15.9, it is bright enough to have been detected by several 
optical and near-infrared sky surveys sinc e 1950 (e.g., 
IZacharias et~alll2004l : ISkrutskie et a"l]l2006f l. iBond et all 
()2002l ) first drew attention to this source following its de- 
tection in the Faint Images o f the Radio Sky at Twenty 
Centimeters (FIRST) survey ()Becker et al.lll995[ ). Their 
follow-up optical observations revealed flickering and a 
spectrum showing strong, double-peaked hydrogen and 
helium emission lines. Based on these, they suggested 
that the system was a cataclysmic variable (CV), that 
is, a binary system in which a white-dwarf primary ac- 
cretes from a close companion through Roche-lobe over- 
flow. Aro und the same time, o ptical spectra obtained by 
SDSS an dlSzkodv et al.l (I2003D appeared similar to those 
found bv lBond et al.l (|2002[ ). 

It was therefore surp rising w hen Thorstensen & Arm- 
strong (2005; hereafter |TAO^ obtained further spectra, 
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starting in 2003 January, that showed substantially lower 
flux and an ordinary-looking G-type spectrum, w i th no 
detectable emission lines (see also iHomer et al.l l2006l 
who found t h e sam e G-type spectrum in 2002 May). 
IWoudt et ai] (|2004l ) published time -series photometry 
that showed a smo oth, ~ 0.3 mag modulation at a 4.75- 
hour period: [TA05| confirmed this result, and also found 
that the radial velocity of the G-type spectrum varies on 
the same period with an amplitude of 268 ± 4 km s^^. 
It has thus been established that J1023 is a binary and 
the G-type spectrum is from the companion star in the 
binary. Furthermore, on the basis o f the r adial velocity 
and results from light-curve fitting, ITA05I inferred that 
the compact primary in J1023 was likely to be a neu- 
tron star or possibly a black hole, and suggested that 
the system resembles a Low Mass X-ray Binary (LMXB) 
rather than a CV. X-ray emission frorn the sy stem was 
subsequently detected by iHomer et al.l ()2006f ). but the 
X-ray luminosity is low, Lx ~ 10^^ ergs s~^ for source 
distance D = 1.3 kpc. This is considerably short of the 
^ 2 L0 compact o bject lu minosity needed to account for 
the heating effect ()TA05| ). 

The puzzles set forth by TA05 were resolved when 
Archibald et al. (2009) discovered pulsed radio emission 
from J 1023 in an untargeted radio pulsar survey carried 
out using the Green Bank Telescope in mid-2007, thereby 
confirming that J1023 is a neutron star binary. More- 
over, the pulsar's spin period P is only 1.69 ms, making 
it the fifth fastest among ~170 k nown milli second pulsars 
(MSP s; P < 30 ms; e.g., Manchester et al.. ,2005: Lorimei] 
l2Ci08h . 

The discovery of the millisecond pulsar makes the 
bright spectrum observed during 2000-2001 extremely 
interesting. The double-peaked emission lines are com- 
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monly seen in CVs and LMXBs, and are a typical feature 
of accretion disks. This suggests that during the period 
in which emission-hne spectra were observed, the com- 
panion star overflowed its Roche lobe and mass transfer 
occurred, forming a disk surrounding the MSP. It is not 
clear whether significant mass accreted onto the neutron 
star, because no evidence of strong, accretion-powered 
X-ray emission has been found. In any case, by showing 
both an accretion disk and radio pulsations, this system 
is likely the first such binary found at the end of its tran- 
sition from an LMXB to a radio MSP. It is a valuable 
laboratory for studying the evolution of X-ray binaries 
and the formation of radio MSPs. 

In this paper, we report the results of the first detailed 
analysis of the SDSS spectrum. We summarize in § 11.11 
the binary properties of J1023, estima ted fro m previous 
optical and recent radio observations (|TA05l : Archibald 
et al. 2009). In §[2] we briefiy describe the SDSS spec- 
troscopic observation of J1023, and present our study of 
spectral features in §[3l From this, we derive the proper- 
ties of the short-lived disk. We discuss the implications 
of the results in § ID 
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Fig. 1. — Mass and radius values (solid curve) derived for 
the companion star in J1023 when Afns = 1-3 Mq is consid- 
ered. Several i values are marked on the curve (diamonds). For 
Mns = 1.4 Mq, M2 = O.2M0 and i = 46°. The m ass-radius rela- 
tion for low-mass main-sequence stars (Tout^et^l.| 1996) is plotted 
(dashed curve), indicating that main-sequence stars cannot fill in 
the Roche lobe. For illustration, we also show the radius values for 
the companion star when the radius is only 85% that of the Roche 
lobe (dotted curve). 

1.1. Binary Properties 

From radio timing of the pulsar in J1023, the pulsar's 
Keplerian mass function is found to be 1.1x10"'^ Mq. 
Combining this with the measured radial velocity ampli- 
tude of the companion star, the mass ratio q is deter- 
mined, q = Mns/M2 = 7.1 ± 0.1. Here Mns and M2 are, 
respectively, the neutron star and companion masses. 
For a canonical neutron star mass Mns = 1-4 Mq, this 
implies M2 — 0.2 Mq. If we allow a wider range of 
Mns from 1 to 3 Mq, then M2 = 0.14-0.42 Mq. Ac- 
cordingly, the orbital inclination angle i is limited to 
34°-53°, nearly identical to the range of inclinations es- 
timated from the heating effect bv .TAOS, . If in addition 
the radius of the companion star R2 is equal to that of 



its Roche lobe^, then the companion star's radius is al- 
most entirely a function of its mass. Figure [1] shows the 
deri ved ma ss and radius values for the companion star. 
The I^Ol light curve fits showed the companion to have 
an effective temperature Toff of 5600-5700 K, similar to 
that of a mid G-type dwarf. Based on the discussion 
above, we adopt canonical values of Mns = 1.4 Mq and 
M2 — 0.2 Mq throughout this paper, unless mentioned 
otherwise. 

2. SDSS SPECTROSCOPY 

The spectrum of J1023 was taken on 2001 Febru- 
ary 1 using the dedicated SDSS 2.5-m te lescope at 
Apache Point Observatory (jYork et al.l l2000f), and was 
included in the 7th SDSS data release (|Abazaiian et al.l 
I2009D . The spectrum covers from 3800 to 9200 
A with spectral resolution A/AA ~2000; the total 
on-source time was over 1 hr. It was reduced us- 
ing the SDSS Spectro2d pipeline (version v5_3_12; 
details can be found on the SDSS webpages, e.g., 
http: / /www. sdss.org/ dr7/products /spectra/index. html) . 
In the pipeline fiux calibration procedure, a few standard 
stars are observed on each spectroscopic plate, and flux 
calibration is achieved by comparing the spectra of the 
standard stars to their model spectra. Generally, the 
3a uncertainty of SDSS spectra in r' filter is 0.15 mag. 
The SDSS spectrum of J1023 is shown in Figure [21 The 
wavelengths are vacuum heliocentric. 

3. SPECTRAL ANALYSIS 
3.1. Emission Lines 

The SDSS spectrum shows emission lines, mostly from 
the hydrogen Balmer series and He I, on a smooth, blue 
continuum. All the emission lines are double-peaked. 
The spectrum resembles those of compact LMXBs and 
CVs; in those systems, the optical emission arises mostly 
in an accretion disk, and the double-peaked profiles re- 
flect the line- of-sight components o f the disk rotation 
velocity (e.g., iHorne fc MarshI Il986f ). The Bowcn flu- 
orescence blend near A4640, which is commonly seen in 
LMXBs, is notably absent. This emi ssion feature consists 
of N III, C III, and O II lin es (e.g.. iKallman fc McCravl 
119801 : [Schachter et al .Hl989t) and is thoug ht to arise froni 
X-ray irradiated c ompanion stars (e.g., iDeguchil 119851 : 
ISteeghs fc Casares|[2002i ). The absence of the A4640 fea- 
ture suggests the lack of X-ray irradiated, highly ionized 
gas in the system and probably indicates a low X-ray 
luminosity at the time. This is consistent with expec- 
tations given the X-ray flux upper limit of the source 
[^2-10 kev ^ 1-2 X 10~^° ergs s~^ cm~^, derived from 
the Rossi X-ray Timing Explorer (RXTE) All-Sky Mon- 
itor data; Archibald et al. 2009]. LMXBs that show the 
Bowen emission feature typically have large X-ray lumi- 
nosities, ~10^^ erg s~^. 

For each emission line, we measured the total flux 
and the wavelengths of the red and blue peaks in 
the following manner. We first dereddened the spec- 
trum using the wavelength-dependent extinction law 
from iFitzpatricS (|1999D and taking Ay ~ 0.14 from 

From studies of MSP binaries, the companion stars h ave of- 
ten been found to be close to filling; th eir Roche lobes (ITA05I : 
IStappers eral|[200ll : IRevnolds e"^aI|[%07^ . 
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Fig. 2. — SDSS spectrum of J1023. The emission lines in the spectrum are la belled. The arrow in t he top panel indicates the absence of 
the A 4640 A emission feature. This feature is commonly seen in LMXBs (e.g.. Ivan Paradiis fc McClintockiil995i ^. in which it arises from 
X-ray irradiated companion stars. 



ISchlegel. Finkbeiner. fc Davis! ()1998f ). We then esti- 
mated the continuum with a 6-order polynomiaL To fit 
the continuum-subtracted double-peaked fines, we used 
either two-Gaussian or two-Lorentz functions. For most 
of the fines, the two-Lorentz profile gave the better fit. 
Tables [1] and [2] give the results for the hydrogen and he- 
lium lines, respectively. The vacuum wavelength values 
used are from the atomic spectra database at the Na- 
tional Institute of Standards and Technology (NIST).* 
Because of the crowding of the lines in the blue region, 
it is difficult to determine the continuum flux exactly; in 
measuring those lines, we estimated the continuum from 
the polynomial fit. 

We focus here on the hydrogen lines, since they are 
stronger than the helium lines and thus more accurately 
measurable. As can be seen, the blue and red peaks of 
the 11/3, 7, e, and C lines have velocities in the ranges of 
-(610-710) and (470-490) km (Figure^. The aver- 
age values are —660 km s^^ and 480 km s~^. Comparing 
to these four lines, US has a stronger, more red-shifted 
profile, while the blue and red peak velocities of Ha are 
significantly smaller. The blue and red peaks of the lines 
are not symmetric, suggesting contributions from other 
components. A hot spot on a disk, arising from the inter- 
action between the gas flow from a compan ion star and 
the disk, can cause suc h asymmetry (e.g., lOrosz et al.l 
[T99llMason et al.l[2000h . 

The fluxes of the Balmer lines are comparable with 
each other, indi cating that the disk was optically thick 
(| Williams! llQSOf) . For example, the peak intensity (in 
frequency units) ratio of Ha to H/3 is ~1.7. This value is 
lower than that found in the classical case of an optically 

® see |http:/ /physics. nist.gov/PhysRefData/ASD/lines-form. html] 



thin gas fe.g. . IOsterbrock!ll974f ) . In addition, the central- 
valley regions of the Balmer lines clearly have a shape of 
'V, not 'U' (FiguresElandlll ) , also indicating an optically 
thick disk (jHorne fc Marsh[|1986l ). The double-peaked 
shape of the lines is generall y considered to be a function 
of sev eral disk parameters |Smak!lT981l : iHorne fc Marsh! 
!l986f ). The line intensity, F{u), can be written as 



F{u) oc 



y,3/2-Q^^ 



(1 



y2j,^l/2 



[l-f4T2(i)w2^(l-M2^)]l/2, 



where r is the normalized disk radius (r = 1 corresponds 
to the outer edge of the disk), u is the dimensionless ra- 
dial velocity (u = 1 at r = 1), ri is the ratio of the inner 
disk radius r-m to the oute r disk radius rput, rz=min(l, 
M~^) and T = sinitani (jSmak! !l98l"l : !Horne fc Marsh! 
!l986l : !Orosz et al.! !l994'). In the equation, a power-law 
model, /(r) oc r~", is assumed for the disk density func- 
tion of the emitting atoms, where a is usually found to 
have values between 1.0-2.0. 

We used this profile function to fit the Balmer lines in 
the SDSS spectrum. The free parameters are peak veloc- 
ity Vp, i, ri, and a. Because the lines are obviously not 
identical and there are relatively large flux uncertainties 
in the blue region, we flt the Ha line mainly to constrain 
the parameters. Although the two Ha peaks are not sym- 
metric, the minimum value is 113 for 108 degrees of 
freedom, indicating that the disk profile fits the Ha line 
reasonably well. The fitting is sensitive to a: a ~ 1.6 al- 
ways provides the best fit. Vp'is tightly constrained in the 



narrow range of 395(— 20, -1-10) km s~ (99% confidence). 
Because changes of ri only cause small disk area and flux 
changes and i is limited by the central valley region, only 
a small fraction of the whole line profile, ri and i can- 
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Fig. 3. — Balmer H/3-Hf emission lines in the SDSS spectrum. The blue- and red-shifted peaks are not symmetric. The right wing of 
the red-shifted HiS peak may be contaminated by another line component. 



not be tightly constrained. The 99% confidence ranges 
are ri < 0.04 and i < 44° (note that since the inner 
disk radius cannot be smaller than that of the neutron 
star, ri has a lower limit, ri > 10~^). We tested fitting 
H/3 and found different results: Vp = 540 ± 30 km s~^, 
n = 0.0002-0.08, and i < 60° (99% confidence; a ~1.6). 
The deep valley in the H/3 line and ot her bluer hydro- 
gen li nes favors a large inclination angle (|Horne fc MarshI 
|1986[ ). However, inclination angles determined from fit- 
ting double-peaked lines obviously are not reliable. For 
example, the hydrogen lines in the SDSS spectrum have 
different central depths, which would suggest different 
inclination angles. Also, in other cases, the central val- 
ley o f a line from the sa me source was seen to be variable 
(e.g.. lOrosz et al.l |1994D. 
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Fig. 4. — Ha emission line in the SDSS spectrum. The disk 
model profile with a = 1.6, Vp = 395 km s~^, ri = 0.004, and 
i = 40°, is plotted as the dashed curve. Because of the asymmetry 
in the two peaks, the model profile does not fit the peak regions 
well. 

If we assume the disk is Keplerian, the disk outer edge 
rotation velocity Vd and rout can be inferred to be Vd = 
Vp/sini and w = (GAWVj) ~ l.l-1.3x 10" sin2(i) 
cm. The Roche lobe radius Ri of the neutron star 
can be found and the disk is expected to be cut off at 



the tidal radius i?tidcs = 0.9i?i = 5.7 x 10^° cm (for 
M2 = 0.2 Mq and M^s = 1-4 Mq). Since Tout cannot be 
larger than i?tides5 we find i < 41°-46°, consistent with 
the low inclination angle value set by the mass ratio. On 
the other hand, this implies that rout — -Rtidcs, because 
when Mns = 1.4 AIq is assumed, i should not be sig- 
nificantly different from 46°. We checked the case when 
Mns = 2.9 Mq. The resulting i is between 31°-34°, also 
approximately consistent with the value given above in 
§ 11.11 The larger peak velocities of other hydrogen lines 
probably suggest that they were emitted from a smaller 
disk area. For example, for Vp = 540 km (H/3 peak 
velocity) and i = 46°, rout ~ 3.3 x 10^" cm. 

The peak velocities of the hydrogen lines have negative 
mean values, which would suggest that the disk (and the 
neutron star) was moving towards us. The velocities of 
the lines are significantly different. If we consider only 
H7, e, and C lines, which have similar velocity values, the 
average radial velocity is —72 ± 36 km s~^, not w ell de- 
termined. This velocity would be orbital because ItA05I 
found nearly zero radial velocity for the binary system. 
Considering the radio timing measurements of the system 
and the orbital phase 0.65-0.92 (phase 0.0 corresponds 
to the ascending node of the pulsar) for the SDSS ob- 
servation, the radial velocity of the neutron star was in 
the range of (— 22)-(-f 33) km s^^ during the observation. 
The median velocity value is ~10 km s~^. These possi- 
ble radial velocity values are different from the observed 
value, suggesting that either the lines were distorted by 
emission from other components such as a hot spot, or 
the disk was not axis-symmetric (e.g.. lMason et aLll2000l 
and references therein) . Generally, radial velocity curves 
derived from emission lines from a disk may not be a 
reliable indicator for the primary's o rbital motion (e.g., 
lOrosz et aLlll994l: iMason et al.ll2000D . 

The helium lines in the SDSS spectrum are rela- 
tively strong. For example, the flux ratio between 
Ha and He I A6680 is approximately 3.3, lower 
than typical va lues found in normal CV systems (e.g., 
iThorstensen et al. 2002 lj). It is like ly that helium is en- 
hanced (j Williams fc Fergusonl[l98^ . As can be seen in 
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Figure [U the companion must be much less massive than 
a main-sequence star of the same radius. Also, it is 
much hotter than a main-sequence star in the allowed 
mass range. While J1023 is not a CV, it is similar in 
that it is a semi-detached binary; as a general rule, sec- 
ondary stars in CVs tend to be slightly cooler th an they 
'shou ld' be based on main-sequence expectations (|Kniggel 
l2006f ). The secondary in J1023 is therefore highly anoma- 
lous. A few CVs do show sim ilar 'too-hot' secondary 
stars (iThorstensen et aD '2002c','a'): in those cases, the sec- 
ondary has evidently undergone significant nuclear evo- 
lution before losing large amounts of mass, leaving an 
undermassive core enriched in helium. The secondary in 
J1023 appears to be a similarly evolved, stripped core, 
much hotter, and much larger in radius, than one would 
expect given its mass. The star p robably had a mass of 
> IMq at onset of mass transfer ()Delovell2008() . and has 
lost much of its mass during its X-ray binary phase. 

3.2. Continuum 

The SDSS spectrum is substantially brighter than 
those spectra obtained after 2002 May; to sho w this, its 
continuum is plotted with the spectrum m 
Figure O The latter spectrum, also dereddened with 
Av = 0.14, is a mean of 23 spectra obtained in 2003- 
2004. Based on the spectra, a mid-G spectral type 
has been identified for the low-mass companion star in 
J1023. The average V magnitude of the companion 
is 17.5 (jTAOSj), corresponding to a dereddened flux of 
4.3x10"^^ ergs~^ s~^ cm~^ A~^. In Figure [H we plot 
a G5V spectrum smoothed and normalized to the flux 
value at V wavelength; the spect rum is from th e stellar 
spectral flux library published bv lPicklei ()1998[ ). As can 
be seen, the library spectrum is generally consistent with 
that of the quiescent emission from J1023 — the blue-end 
region of the former is lower than that of the latter, but 
the flux uncertainties of the latter in the region are as 
large as 15%. A more stringent test of how well the 
spectrum matches a G5 dwarf could be made if we had a 
flux-calibrated spectrum of J1023 covering wavelengths 
to 3500 A, as a sharp turnover is expected at 4500 A. In 
the red end, if the G5V spectrum were extended to near- 
infrared JHKs wavelengths, it would have flux values 
approximately consistent with the 2MASS measurements 
of J1023 (J = 16.3±0.1, H = 15.7±0.1, Ks = 15.9; the 
uncertainty on Kg could not be determined, and the ob- 
servations were carried out on 2000 February 6). All in 
all, the evidence supports the use of the G5V spectrum 
for the companion star. 

Including the normalized G5V spectrum, we fit the 
SDSS continuum with a simple disk model. The fiux 
at wavelength A from the disk is 

A-kHc^ cos i /■''°"' rdr 

where h is the Planck constant, c is the speed of light, 
and k is the Boltzmann constant. We assumed that the 
disk temperature Td is a power-law function of the disk 
radius r, Td oc r~P , and further considered p = 3/4, the 
same as that in the steady thin disk model ([Frank et al.l 
Il992l ). If the neutron star had strong X-ray emission at 
the time due to accretion, which is unlikely based on the 
X-ray upper limit (Archibald et al. 2009) and absence 



of the B owen fluorescence f eature, p = 3/7 might be 
assumed (jVrtilek et al.lll990f ). Under the simple model, 
the more likely case would be that p varied during the 
disk lifetime. If that is the case, multiple observations of 
a future disk phase would be needed to detect spectral 
changes caused by disk evolution. 

We first tried three parameters in the fitting: D, rin, 
and the temperature at rjn, while Tout was set to be 
rout = 0.9i?i « 5.7 X IQi" cm (for M„s = 1.4 M© and 
i = 46°). The distance D is constrained by the flux 
from the companion, since it is a function of (i?2/£')^- 
For example, if we assume a G5V dwarf (Toff ~ 5600 K, 
radius i?* = 0.92 Rq, and absolute magnitude My = 
5.1) as the companion, V — 17.5 would imply a distance 
of 2.8 kpc. When i?2 is assumed to be equal to the Roche- 
lobe radius (0.38 Rq) of the companion, D ~ 1.2 kpc, 
inferre d from {R2ID) = (O.92i?0/2.8 kpc). Similarly, 
lTA05l have found D ~ (2.2 kpc)(M2/Mo)i/3, and D ~ 
1.3 kpc when M2 = 0.2Mq. Therefore, we searched for 
the best fit in the range oi D = 0.9-1.7 kpc. 
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Fig. 5. — Comparison of the SDSS continuum from 2001 (top 
liglit dashed liistogram over-plotted with fiux uncertainties short 
lines) with the mean J1023 spectrum (bottom solid curve) obtained 
in 2003-2004. The latter has substantially lower flux, and can 
be described by a G5V spectrum (dash-dotted curve). The error 
bar at the left side indicates the large uncertainties in the blue- 
end region of the mean spectrum. Combining the G5V spectrum 
with that of a disk model (dotted curve), the summed spectrum 
(top dark solid curve) fits the SDSS spectrum well, supporting the 
existence of a disk during 2000-2001. 

Optical emission from J1023 is modulated because the 
visible area of the heated face of the co mpanion star 
varies as a function of orbital phase (|TA05l ) . We assumed 
this flux modulation the same in the bright state as in 
the quiescent, which is expected since no strong X-ray 
emission (to heat the companion) was seen. The SDSS 
observation was made at orbital phase 0.65-0.92, right 
within the fl at top of the companion's modulated light 
curve ()TA05l : note t hat bec ause phase 0.0 is defined dif- 
ferently here from in lTA05l . the phase in this paper leads 
by 1/4). The average V magnitude during the phase 
range is 17.354±0.014, implying a 15% flux increase at 
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the time from the companion. We simply increased the 
flux of the G5V spectrum by 15%. Because the flux 
increases are caused by the heated half surface of the 
companion star, they are actually wavelength-dependent 
(TA05). However, since the SDSS spectrum is much 
brighter than that of the companion, the differences are 
negligible in the fitting. 

The flux uncertainties on the SDSS spectrum are ap- 
proximately 2% in the middle region and 4% at the end 
regions. If we use the uncertainties, the minimum re- 
duced (2922 de grees of freedom) resulting from the 
fitting is 4.6. We therefore included a systematic uncer- 
tainty, which on average could be ~5% in the r' band. In 
addition, it should also contain uncertainties on the com- 
panion's spectrum, since no uncertainties for the G5V 
library spectrum were assumed. This systematic uncer- 
tainty was added in quadrature, while its value was ad- 
justed to have the minimum reduced = 1. A value 
of 3.4% for the uncertainty was required. The fitting is 
sensitive to the parameters, and we found D ~ 1.0 kpc, 
Hn ~ 1.5 X 10^ cm, and ~ 34000 K. The ri„ value 
is consistent with those from our Ha fitting, while D is 
~20% lower than those derived above and bv lTA05l . The 
D value may suggest that i?2 was not equal to, but lower 
than the Roche lobe radius, li D = 1.3 kpc is required, 
the spectrum is less well fit, and rjn ~ 1.4 x 10^ cm and 

— 40000 K. In Figure El the best-fit model spectrum 
is shown. 

We further considered that this was a steady thin 
disk case. Given the values obtained above, the mass 
accretion rate M in th e disk can be est imated from 
Td ~ 8000Mi/Vj"o^/'' K (iFrank et al.lll992D . Mig ~ 1.1, 
where rig and Mig are r and M in units of 10^" cm 
and lO'^^ g s~^, respectively. Based on this M value, 
the Alfven radius tm, which is expected to be the in- 
ner disk radius, is much smaller than the obtained 
ru 2.5 X 106(a'/i6/1.1)-2/7(_b/io8 Qf^ cm, where 
B is the surface magnetic field strength of the neutron 
star (B ~ 10^ G; see §2] below). The discrepancy indi- 
cates inconsistency in the results when the steady thin 
disk case is considered. We therefore tested to find the 
alternative best fit by setting Tin ^ Rvi — 2.4 X 10 cm 
(where i?n is the corotation radius) and Td oc M^/'^r~^/^ . 
Generally, when ~ i?f2, accretion onto a neutron star 
occurs and the pulsed radio emission from the pulsar 
is quenched due to accretion. It can be noted that to 
satisfy Tin ~ tm ~ i?n, Miq ^ 1. Since tm is weakly 
dependent on M, a large range of M values will result 
in small changes of tm (or rin), indicating that it is rea- 
sonable to fix Tin at the Rq, value. To search for a good 
fit, Tout was set to be a free parameter in the reasonable 
range: rout < 5.7 x 10^" cm. We found from the fitting, 
M ~ 1.0 X 10^^ g s~^ (which is consistent with the as- 
sumed Tin value), Tout = 5.7 x 10^° cm, and D = 1.1 kpc, 
but the minimum reduced — 1-9- In order to have 

— 1, a systematic uncertainty of 5% would have to be 
included in the fitting. 

4. DISCUSSION 

We analyzed the SDSS spectrum, which was obtained 
in 2001 February when J1023 was in a bright state. Al- 
though the double-peaked hydrogen lines in the spectrum 
have shapes more complex than those in a standard disk 



profile, our detailed study of them shows that their prop- 
erties are consistent with those of a disk around the neu- 
tron star in J1023. Considering a G5V spectrum normal- 
ized to the brightness of the companion star, we have also 
found that the SDSS continuum can be fit with a sim- 
ple disk model, supporting its disk origin. These studies 
thus help indicate that in the bright state a disk existed 
in the binary, demonstrating the important feature im- 
plied by the source: at the beginning of a radio MSP life, 
its companion star in the binary is still able to overflow 
its Roche lobe and a disk can thus be formed around the 
MSP. 

From fitting the continuum, the temperature profile of 
the disk is estimated to be 9100r^Q'^^^ K. Such a tem- 
perature profile implies M ~ 1.1 x 10^^ g s^^ and a 
small, ~ 2.5 x 10^ cm Alfven radius when the standard 
steady thin disk is assumed. The radius value is much 
smaller than the rjn value obtained from the fitting, in- 
dicating inconsistency in the results. However, it can 
be shown that given the derived M value, the viscous 
(or radial drift) timescale could be as long a s ^110 days 
(the viscosity coefficient a = 0.1 is assumed: iFrank et al.l 
Il992f) for the outer disk. The starting time of the disk 
formation is constrained by the SDSS imaging observa- 
tions, which were made on 1999 March 22 and indicate 
that J1023 was in the quie scent state at t he time (e.g., 
q' = 17.99, r' = 17.43; Bo nd et al.l[200l : ISzkodvetaLl 
120031) . The disk formed at some time later but before 
2000 Ma y 6 when the first bright-state spectrum was 
taken by iBond et al.l (120021 ) , and had existed for > 8 
month but < 2 yr before the SDSS spectroscopic ob- 
servation. Because the viscous timescale is comparably 
long, these constraints suggest that the disk could have 
had Tin ~ 10^ cm if it had not fully extended to the 
allowed smallest radius. We note that B currently is es- 
timated to be low, B ~10^ G (although the spin-down 
rate P of the pulsar is quit e uncertain, P ~ 7 x 10"^^ 
s s~^: [Archibald et al.ll2009[ ). supporting a small Alfven 
radius. 

Our assumption of a standard thin disk could be too 
simplified. On the other hand, it is plausible that the 
SDSS spectrum was not well calibrated in fiux — the real 
spectrum would be steeper and thus rjn could be as small 
as ~ 2.5 X 10*^ cm (Mie = 1.0 and B = l{f G are used). 
This Tin value would be comparable to i?si, suggesting 
that accretion onto the pulsar could have occurred. How- 
ever, the M value would imply an accretion luminosity 
of Lx ^ GMnsM/Rns =i 2 X lO^s ergs s^^. Comparing 
this luminosity to the X-ray flux upper limit, either the 
source is >10 kpc away or no accretion onto the neutron 
star occurred. Since multiple lines of evidence point at a 
distance of ~1.3 kpc, the latter is the likely case. 

The fight cylinder (LC) radius of the MSP in J1023 
is He ~ 8.1 X 10^ cm. Generally, if Rq 4C Hn < 
He, t he system would ha,ye be en in the propeller 
phase (jlllarionov fc Sunvaevlll975( ) and the radio pulsar 
would have been que nched by the accretion flow (e.g., 
ICampana et~anil998l) . The minimum luminosity Lie pro- 
duced from the accretion flow would have been when 
Hn ~ nc, Lie ^ GMnsM/ric ~ 2 X 1035Mi6 erg s'^. For 
D = 1.3 kpc, Pmin ^ lO-^Mie erg s'^ cm^^. This flux 
would have been detected by the RXTE All-Sky Moni- 
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tor. The upper limit on M given by the X-ray flux upper 
hmit is 10^^ g s~^. Therefore, if the M value we have 
derived is approximately correct, rjn would have been 
larger than r\c and the radio pulsar would not have been 
quenched at the time. This is consistent with the large 
rin value we obtained. On the other hand, if the M value 
is overestimated, the disk could have extended inside the 
LC and quenched the radio pulsar. 

We estimate the total disk mass Md by simply us- 
ing the standard, a-prescription disk model and in- 
tegrating the surface density E over the disk area. 
In the model, S ~ h.la-'^/^Ml|^^r~^''^ g cm^^ - 
36(a/0.1)~^/'''Mjg"'^°rj~Q'^^* g cm"^. Using the inner and 
outer disk radii obtained above, Md ~ 1.7 x lO^'^ g. The 
formation time would have been 0.5 yr if the mass trans- 
fer rate had been constant and equal to 10^^ g s^^. This 
timescale is consistent with the constraints on the disk 
formation time set by the optical observations. 

The existence of a disk for more than a year makes 
J 1023 an interesting case for studying the pulsar-disk 
interactio n. According to the s tandard pulsar accretion 
scenario (jCampana et al.lll998l and references therein), 
the minimum mass accretion rate M-^ai-a in a disk, which 
is required for the disk to be stable against a pulsar's 
radiation pressure, is obtained when tm ~ r\c. Detailed 
calculations made by lEk^I fc Alparl ()2005D show that a 
disk can still be stable if its inner radius is only slightly 
larger than n^- For the J1023 MSP, Afmin 2 x lO^"^ 
g s""^ . With such a low rate, the disk in J1023 would 
have lost most (~95%) of its mass before the final disk 
disruption. The mass loss rate of the disk would have 
been 3x10^^ g s~^ on average, where 2.0 yr is assumed 
for the disk lifetime. Because the spin-down luminosity 
of the pulsar is probably as l arge a s ~ 5 x 10'^^ erg s~^ 
(from the light curve fitting, iTAOSi have found that the 
required luminosity from the MSP is ^ 10^'' erg s~^, 
supporting the spin-down luminosity value), the mass 
loss could have been due to the pulsar wind ablation of 
the disk (|Miller fc HamiltonI 120011 ). On the other hand, 
if the disk had extended inside the LC and quenched 
the pulsar wind, the mass loss would have been caused 
by the propeller effect: the mass inflow is halted by the 
magnetic field outside of the corotation radius and the 
material is repelled away from the system. 

It is likely that J1023 will repeat its process of having 
an outflow and forming a disk, although the timescale be- 
tween two such events is not known. If we do see the next 
event in the near future, it would provide a great oppor- 



Abazajian, K., et al. 2009, ApJS, 182, 543 

Archibald, A. M., et al. 2009, Science, 324, 1411 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 

Bond, H. E., White, R. L., Becker, R. H., & O'Brien, M. S. 2002, 

PASP, 114, 1359 
Campana, S., Colpi, M., Mereghetti, S., Stella, L., & Tavani, M. 

1998, A&A Rev., 8, 279 
Deguchi, S. 1985, ApJ, 291, 492 

Deloye, C. J. 2008, in American Institute of Physics Conference 
Series, Vol. 983, 40 Years of Pulsars: Millisecond Pulsars, 
Magnetars and More, ed. C. Bassa, Z. Wang, A. Gumming, & 
V. M. Kaspi, 501-509 

Ek§i, K. Y. & Alpar, M. A. 2005, ApJ, 620, 390 

Fitzpatrick, E. L. 1999, PASP, 111, 63 



tunity for studying the pulsar-disk interaction. The short 
lifetime and brightness of the disk would allow close mul- 
tiwavelength monitoring, from which we might be able 
to test current pulsar accretion theory. For example, we 
might find at which stage of the disk evolution pulsed 
radio emission is quenched. Would it be, as generally 
suggested, when the disk extends right inside the LC? If 
the transition to the propeller phase occurs, we might be 
able to seek direct evidence for mass loss due to the pro- 
peller effect (e.g.. iRomanova et al.1120 09'). On the other 
hand, we would study the evolution of an accretion disk 
as its being ablated by a pulsar wind. More importantly 
we might directly observe how a disk is disrupted by a 
pulsar. 

Funding for SDSS and SDSS-II has been provided by 
the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. 
Department of Energy, the National Aeronautics and 
Space Administration, the Japanese Monbukagakusho, 
the Max Planck Society, and the Higher Education 
Funding Council for England. The SDSS Web Site is 
http: / / www.sd ss.org/[ 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Par- 
ticipating Institutions are the American Museum of Nat- 
ural History, Astrophysical Institute Potsdam, Univer- 
sity of Basel, University of Cambridge, Case Western 
Reserve University, University of Chicago, Drexel Uni- 
versity, Fermilab, the Institute for Advanced Study, the 
Japan Participation Group, Johns Hopkins University, 
the Joint Institute for Nuclear Astrophysics, the Kavli 
Institute for Particle Astrophysics and Cosmology, the 
Korean Scientist Group, the Chinese Academy of Sci- 
ences (LAMOST), Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck- Institute for Astrophysics (MPA), New Mexico 
State University, Ohio State University, University of 
Pittsburgh, University of Portsmouth, Princeton Uni- 
versity, the United States Naval Observatory, and the 
University of Washington. 

This research was supported by NSERC via Discovery 
Grants and by the FQRNT and CIFAR. VMK holds a 
Canada Research Chair and the Lorne Trottier Chair in 
Astrophysics & Cosmology, and is a Fellow of the Royal 
Society of Canada. JRT acknowledges support from the 
U.S. National Science Foundation, through grants AST- 
0307413 and AST-0708810. IHS acknowledges support 
from the Swinburne University of Technology Visiting 
Distinguished Researcher Scheme. 



Frank, J., King, A., & Raine, D. 1992, Accretion Power in 
Astrophysics. (Camb. Astrophys. Ser., Vol. 21,) 

Homer, L., Szkody, P., Chen, B., Henden, A., Schmidt, G., 
Anderson, S. F., Silvestri, N. M., & Brinkmann, J. 2006, AJ, 
131, 562 

Home, K. & Marsh, T. R. 1986, MNRAS, 218, 761 
Illarionov, A. F. & Sunyaev, R. A. 1975, A&A, 39, 185 
Kallman, T. & McCray, R. 1980, ApJ, 242, 615 
Knigge, C. 2006, MNRAS, 373, 484 

Lorimer, D. R. 2008, Living Reviews in Relativity, 11, 8 
Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 
129, 1993 

Mason, E., Skidmore, W., Howell, S. B., Giardi, D. R., Littlefair, 

S., & DhiUon, V. S. 2000, MNRAS, 318, 440 
Miller, M. G. & Hamilton, D. P. 2001, ApJ, 550, 863 



8 Wang et al. 



TABLE 1 
Hydrogen emission lines in J1023 



Line 


Avac 


5X 


V 


FWHM 


EW 


Flux/10-18 




(A) 


(A) 


(km s-i) 


(km s-i) 


(A) 


(ergs cm"^) 


HC b 


3890.1 


-8.9±0.3 


-690 ± 23 


610±210 


3.01±0.22 


78.0±5.7 


HC r 


3890.1 


6.3±0.4 


488 ± 31 


720±260 


2.76±0.23 


71.1±5.8 


He b 


3971.2 


-9.4d=0.3 


-709 ± 22 


640±210 


3.07±0.21 


77.1±5.3 


He r 


3971.2 


6.3±0.3 


479±24 


660±220 


2.96±0.21 


74.0±5.2 


US b 


4102.9 


-9.1±0.3 


-669 ± 22 


690±210 


2.14±0.13 


51.2±3.0 


m r 


4102.9 


10.0±0.3 


728±25 


1090±320 


3.53±0.16 


83.7±3.8 


H7 b 


4341.7 


-9.3 ±0.2 


-639 ± 16 


600±190 


3.89±0.21 


85.2±4.5 


H7 r 


4341.7 


7.2±0.3 


495±24 


890± 290 


4.76±0.23 


103.4±5.1 


H/3 b 


4862.7 


-9.9 ±0.3 


-610 ± 20 


780±250 


5.09±0.25 


93.0±4.5 


H/3 r 


4862.7 


7.7±0.4 


474±23 


760±250 


4.39±0.23 


79.9±4.2 


Ha b 


6564.7 


-9.5 ±0.4 


-432 ± 19 


770±230 


8.45±0.35 


97.5±4.1 


Ha r 


6564.7 


7.8 ±0.3 


355 ± 15 


600±190 


7.45±0.29 


85.5±3.3 



TABLE 2 
Helium emission lines in J1023 



Line 


Avac 


5\ 


V 


FWHM 


EW 


Flux/lO"!'^ 




(A) 


(A) 


(km s-l) 


(km s-i) 


(A) 


(ergs s-i cm"^) 


Hel 6 


4027.3 


-6.3±0.5 


-472 ± 34 


910±300 


3.22±0.22 


79.2±5.5 


Hel r 


4027.3 


10.1±0.8 


755 ± 58 


580±300 


0.96±0.18 


23.4±4.4 


Hel 6 


4472.7 


-8.1±0.7 


-541 ± 48 


580±290 


1.33±0.29 


27.8±6.0 


Hel r 


4472.7 


4.9±1.2 


327±83 


1500±530 


5.22±0.33 


108.2±6.9 


He II b 


4687.0 


-10.2±0.3 


-650 ± 21 


430±160 


1.69±0.17 


32.7±3.3 


He II r 


4687.0 


9.2±0.8 


591±49 


1340±470 


3.99±0.26 


76.8±5.0 


He I b 


4923.3 


-9.1 ±0.6 


-550 ± 33 


640±240 


1.94±0.20 


34.9±3.5 


He I r 


4923.3 


8.5±0.5 


519±29 


520±200 


1.74±0.18 


31.0±3.2 


He I b 


5017.1 


-7.6 ±0.6 


-452 ± 35 


700±260 


2.29±0.22 


39.8±3.8 


He I r 


5017.1 


11.5±0.5 


691±29 


870±290 


3.98±0.23 


68.9±4.0 


He I b 


5877.3 


-7.6 ±0.6 


-388 ± 31 


720±240 


3.82±0.26 


52.4±3.6 


He I r 


5877.3 


5.9 ±0.3 


300 ± 17 


350±130 


2.40±0.19 


32.8±2.6 


He I b 


6680.0 


-10.7± 1.0 


-482 ± 45 


800±300 


2.80±0.27 


31.2±3.0 


He I r 


6680.0 


9.5 ±0.9 


427 ± 40 


630±260 


2.25±0.23 


25.0±2.6 


He I b 


7067.1 


-10.0 ± 1.0 


-425 ± 44 


620±250 


1.98±0.25 


20.1±2.6 


Hel r 


7067.1 


9.1 ±0.9 


388 ± 39 


560±240 


1.99±0.23 


20.1±2.4 
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